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Abstract 

The generalization of the BGK relaxation model to the special rel- 
ativity setting is revisited here. We deal with several issues related to 
this relativistic kinetic model which seem to have been overlooked in the 
previous physical literature, including the unique determination of associ- 
ated physical parameters, classical, ultra-relativistic and hydrodynamical 
limits, maximum entropy principles and the analysis of the linearized op- 
erator. 



1 Introduction 

Our aim is to give a mathematical description of a gas in certain relaxation 
regimes. This will be done in terms of the relativistic BGK equation. In that 
framework, such gases are regarded as consisting of many microscopic structure- 
less particles; our description will be given in terms of the relativistic kinetic 
phase density. This object allows us to compute tensorial moments giving the 
local macroscopic physical quantities of the gas (particle density, pressure, etc). 
We describe the construction of the relativistic BGK system, analyze the hyper- 
bolic hydrodynamical limits towards the relativistic Euler equations and study 
the well-posedness of the relativistic linearized BGK system, together with clas- 
sical and ultra-relativistic limits. 

The models of kinetic theory describe the time evolution of a collection of 
particles. The main focus of kinetic theory has been on classical particles and on 
the special relativistic framework, but also there arc mathematical advantages in 
using kinetic models in general relativity instead of their hydrodynamic counter- 
parts, especially from the point of view of the ensuing singularities [21 122 [M]. 
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A few years before Einstein established the principles of general relativity in 
1916, Jiittner [37] gave in 1911 the first step in extending the kinetic theory of 
gases to the relativistic context, proposing a generalization of the Maxwcllian 
distribution function which is widely accepted nowadays. The next step was 
due to Walker [HI] who derived the relativistic evolution equation in the ab- 
sence of molecular collisions in 1935. Lichncrowicz and Marrot [44] provide a 
first complete relativistic generalization of the Boltzmann equation (which was 
proposed in 1872) in 1940. The analysis of the Cauchy problem in relativistic 
kinetic theory was first addressed by Choquet-Bruhat [HJ (32] [23], and sub- 
sequently in a wide literature, see for example [2H1 US]- The analysis of the 
linearized relativistic Boltzmann equation was first made in [55], see also [23] . 
The Cauchy problem for the relativistic Boltzmann equation has been recently 
revisited by Strain |58| ; he shows that an initial datum starting close enough to 
a global equilibrium launches a unique global solution which decays with any 
polynomial rate towards that global equilibrium. 

Some hydrodynamic extensions of relativistic models have been obtained 
either using the Chapman-Enskog method or the moment method of Maxwell 
and Grad [30] , see Israel [35] [36] , De Groot et al [33] , Marie [47] Degond et al 
[3] and the references therein. We highlight here that the viscous (parabolic) 
limits frequently lead to an infinite wavefront speed for the transport processes, 
which bear some inconsistency in a relativistic context, where the propagation 
speed cannot overcome that of light [5D]. The hyperbolic asymptotic limits 
have the advantage of eliminating this spurious infinite velocity of the wave 
front propagation and that makes them in general more suitable than parabolic 
limits in relativistic dynamics. In [57] the local-in-time hydrodynamic limit 
of the relativistic Boltzmann equation is analyzed using a Hilbert expansion 
and following the ideas developed in the classical case [H] [T3]. As a result, 
local solutions to the relativistic Boltzmann equation near the local relativistic 
Maxwcllian arc constructed, using a class of solutions to the relativistic Euler 
equations that show up in the hydrodynamic limit. In this hydrodynamical 
context, one of the goals of our work will be to deduce the hyperbolic asymptotic 
in terms of the different parameters involved in the relativistic BGK equation. 

The BGK (Bhatnagar, Gross and Krook) model [TU], proposed in 1954 
for classical particles, became the most important model to solve the integro- 
differential Boltzmann equation. The non-linear quadratic (binary) collision 
term of the Boltzmann equation is replaced in the BGK model by a seemingly 
simpler term (however the nonlinearity becomes exponential), which makes the 
derivation of the transport equations for macroscopic variables easier. A prob- 
lem which can be naturally addressed using the BGK model is that of the 
relaxation of a state of a fluid to equilibrium. 

Our description -which is essentially based on the Marie model SB] - will 
come in terms of a relativistic extension of the hitherto known as BGK model 

df 

being Mf a local Maxwcllian implicitly defined by the requirement of having the 



2 



same moments as the distribution function /, which depends on time t, position 
x and velocity v. From the physical point of view, the density of particles is 
assumed to converge to an equilibrium represented by a Maxwcllian function 
of the velocity v when the time t becomes large. During the last 20 years, the 
BGK model has also found an important application: the derivation of numerical 
schemes, namely kinetic schemes to solve hyperbolic conservation laws, see [53] . 
Some approaches to the relativistic BGK models have been proposed in the 
literature, see for example [TTJ [H ESJ 21 EE] and the references therein. It is 
also remarkable the applicability of this type of models in gravitation in order 
to analyze dispersion properties or the stability of special configurations such as 
galaxies. In fact, the BGK-equilibria (or BGK waves in ID) have been analyzed 
for plasmas and gravitation as stationary solutions of coupled systems of type 
Vlasov, Vlasov-Boltzmann, Einstein- Vlasov, Landau-Boltzmann, Boltzmann- 
Maxwell, . . . , see H H QU CGfl HI HI HQ 131 EH1 023 an the references therein. 

The contents of the paper are as follows: Physically meaningful quantities 
and local equilibria arc introduced in the first Section. A generalization of the 
classical BGK model is introduced then in Sectiond whose behavior in the non- 
rclativistic limit we analyze in Section [33] Then we elaborate on the conserva- 
tion laws of the model and how do they give rise to the relativistic and ultra- 
relativistic Euler equations in some limiting regimes, being this the content of 
Sections 13.21 and [ 3.31 The remaining Sections 2] and [5] proceed with the analysis 
of the linearization of the relativistic BGK operator, culminating with a global 
existence proof for the linearized BGK model. Several computations concerning 
moments of Jiittner distributions and other technicalities are displayed in an 
Appendix. 

1.1 Some notes about special relativity 

The fundamental physical theory involved in this description that we want to 
set is special relativity. Just before detailing its precise role we will introduce 
the related notations and conventions that we will follow during the exposition. 
Some background on special relativity and relativistic kinetic theory can be 
found for instance in H [H [31 EJ [33] . 

The space-time coordinates in the four-dimensional Minkowsky's space are 
n = 0, 1, 2, 3, with x° = ct for the time and x 1 , x 2 , x 3 for the position; here 
c is the speed of light. The metric tensor g^ v and its inverse g^ v are given by 

9iiv = g^ v = 1 if fi = v = 0, —1 if /U = v = 1, 2, 3 and if v ^ u. 

Our greek indices run from to 3 and our latin (spatial) indices do from 1 
to 3. We will use Einstein's summation convention, meaning that any index 
that appears twice in an expression, is understood to be summed over its whole 
range. With the aid of the metric tensor we can perform the operations of raising 
and lowering indices. That is, for any four-dimensional vector u (four-vector 
hereafter), 

gavU 1 ' = u a and g au u u = u a . 
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This works in the same way for general tensor objects. 

We will also consider vectors in the Euclidean three-dimensional space, which 
we will always denote by bold characters. Then, notations like |u| = Vu^u* and 
u • q = u^q* stand for the euclidean norm and scalar product respectively. 

1.2 Microscopic and macroscopic quantities 

We are going to describe the microscopic state of a relativistic gas by means 
of kinetic theory. For that, let us introduce the relativistic phase density 
f(t,x,q) > 0, which represents the density of particles with given spacetime 
position (t,x) and momentum q <G R 3 . We will consider that all the gas parti- 
cles have the same mass m. Then the energy-momentum four-vector is defined 
as 

^ = ( g °,cq), q° : = <V(mc)2 + |q| 2 - 

Note that has energy dimensions, whereas q has momentum units. 

Next we can define some macroscopic moments and the entropy four- vector. 
All these make sense at a given point x^ as soon as f(t, x, •) > is not identically 
zero. The fact that the proper volume element dq/q° is invariant with respect 
to Lorentz transformations (i.e. isometries of the Minkowsky space) [TO] is a 
key physical feature of these definitions. 

1. Particle-density four- vector 

N»{t,x)= [ <f/(M,q)^, 
Jr 3 <7 

2. Energy-momentum tensor 

m J K 3 q° 

3. Entropy four- vector 

m J R 3 V V J 1° 

being kg the Boltzmann's constant and r] = m/h 3 , with h the Planck 
constant. 

We may use the macroscopic moments N fJ -,T fJ,lJ of the relativistic phase 
density / in order to compute other macroscopic quantities of the gas, namely 

1. The proper particle density n/ 

n f = (1) 
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2. The dimensionless velocity four- vector Uf 

n f u) = N". (2) 

Note that u^uj)^ = 1 and then = (y/l + |u| 2 , u). Thus for /i = 0, we 
can deduce the relation 

n f Jl + |u/| 2 = / /(t,a,q)dq. (3) 

3. The proper energy density e y 

e/ = (uf) fJ ,(u f ) v T f " / . 

4. The proper pressure pf 

Pf = l((uMn f )u-9^)T^. (4) 

5. The proper entropy density aj 

1.3 The Jiittner equilibrium 

The generalization of the classical global Maxwellian to this setting is the so- 
called Jiittner equilibrium (or relativistic Maxwellian) [37l [38]. The Jiittner 
distribution can be regarded as a function J(n, f3, u; q) describing the state of a 
gas in equilibrium, depending on five parameters: n > 0, j3 > and u £ R 3 . It 
is given by the formula 

J(n < ft U; q) = (mc)3M (/ 3) eXP {~^ U ^} 

or equivalently 

J(n ' * U; q) = ( TO c)3M(/3) eXp (Vl + HV(- C ) 2 + |q| 2 - u ■ q) } . 

Since J(n, /?, u; q) is thought of as an equilibrium distribution, then n is inter- 
preted as its particle density, u as the spatial part of the four-velocity u (and 
as such u^u^ = 1) and mc 2 / (ksP) as the equilibrium temperature. Note also 
that /3 is dimensionless as M{0) is, which takes the form 

M{fi) = -^-3 / CXp(-A^( TOC )2 + | q |2j dq 

= f cxp\-(3^TTW}dp- (5) 
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In such a way, we have that 



nu" = / g*J(n,0,u;q)^f. (6) 

As regards the size of the exponent in the above formulae, we point out that 

> mc 2 (7) 

which is an straightforward consequence of the Cauchy-Swartz inequality for 
timelike vectors in Minkowski space. 



2 From relativistic Boltzmann to relativistic BGK 

The definitions from the previous section apply to any relativistic phase density 
/ = f(t,x,q), which, in case that no forces are involved, should obey a kinetic 
equation of the following form 

^ = C(/). (8) 

We have a transport part on the left-hand side and a collision part C(f) on the 
right-hand side, as in the non-relativistic setting. In the simplest case particles 
interact only by means of elastic collisions -no other forces are assumed to play 
a role-. Then C(f) should be determined in such a way that the conservation 
laws for particle number, energy and momentum hold: 

^=0, ^=0. (9) 
dx» dx v w 

The relativistic Boltzmann equation is a paramount example of such kinetic 
models. It arises as a particular case of (|SJ) with C(f) = Q{f, /), which stands 
for a bilinear operator in /, verifying that Q(J, J) = 0. The detailed descrip- 
tion of the relativistic collision operator Q(f,f) is inessential for our present 
purposes; we refer the reader to [51 119| [55] for more information of the subject. 
We just point out that the relativistic extension of the BGK model and related 
descriptions -see |45j for instance- arise as model equations for the relativistic 
Boltzmann description. 

2.1 The BGK Model 

Now we are going to construct a BGK-type model which satisfies the conserva- 
tion laws ©. In this way we will recover a model previously posed by Marie 
[151 HE] . Making a parallel with the relativistic Boltzmann equation and the 
classical BGK model, we want to build an equation for / of the following form 

¥ + C fVJ^(J r /):^, (10) 
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the right hand side being a relaxation operator constructed from a local rela- 
tivists Maxwellian Jf, which we will detail below. Here u> denotes the collision 
frequency and q = q/ \J (mc) 2 + |q| 2 . Note that such a model has also the form 
®, with C(f) = Qbqk{I) = mcu(J f - /). 

We define Jf implicitly as a local Juttner distribution J(n, /3, u; q) whose 
parameters n, u, (3 will be given in the next Lemma Q] through some integrals 
of / in an unique way. In order to do that, we impose that functions verifying 
(fT0|) should also satisfy the conservation laws (j9]). That requirement turns up 
to be (via multiplication by 1 and q^ /m and integration in (|10|)) equivalent to 

/''§-/ f% <"» 

fj '?-L"$- (12) 

Actually, let us see that these five conditions allow us to determine the five 
parameters of the Jiitner equilibrium and to construct the relaxation operator. 

Lemma 1. Let f = f(q) > 0, for q G R 3 , a non-zero phase density (in the 
almost everywhere sense) for which the moments N tJ, ,T tJ,l/ exist. Then, we can 
find a unique set of quantities n, u ; j3 and a function Jf := J(n, /3,u;q) such 
that (jlip and (|12p are fulfilled. Moreover, for this choice of Jf , a relativistic 
phase density f verifying (|10|) fulfills the conservation laws 

Proof. We first note that taking n — n/ and u = u/, P^|) is a direct consequence 
of © and ©. 

It only remains to determine the parameter j3 in order that holds. Fol- 
lowing Lemma [T2l condition (|lll) is equivalent to the equation 

K^^rru^ f dq 

K 2 ([3) n f 7 M 3%°' { ) 

The function on the left hand side is known to be increasing and with range 
[0, 1), see Section IA~4l in the Appendix and pQ. It only remains to prove that the 
right hand side is in the same range. In order to do that we use the invariance 
of nf and -§ under Lorentz transformations (see Lemma [5] in the Appendix) to 
write 



/ 4 

Jr3 q u 



/a dq 



where we have considered the Lorentz boost A such that it/ A = (1,0, 0,0). We 
can conclude by noticing that q° — c^J (mc) 2 + |q| 2 > mc 2 . □ 

For future reference we denote /3/ the unique parameter /3 verifying equation 
(fT3)l . We also remark here (see Lemma |8] in the Appendix) that this quantity is 
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Lorentz invariant. So finally the BGK-type model that we propose is given by 
(TO]) together with 

J f = J(n / ,/3 / ,u / ;q). 

The well-posedness of (|TU)) remains as an open problem (compare with [5T1 IB]); 
see however [52] for the 2-dimensional case. 

2.2 Conservation laws 

Under no external forces and no radiation, any matter model should obey the 
constitutive relations (|9]). For our model, these conservation laws can be ob- 
tained from (|10[) upon multiplication by (1, /to) and integration in M 3 against 
^Jr, by noticing that Lemma[T]allow us to ensure that the corresponding integrals 
of the rhs are zero. Then, the conservation laws read: first, 

[ fdq+cdiv x f ^ = dq = 0, (14) 

ot Jk3 J R 3 yj(mc) 2 + |q| 2 

as the equation for (conservation of the number of particles) which, using 
([2jl and (HJ), leads to 



d_ 

ot 



"j 



l + u 2 j + cdiv x (nfUf) = 0. 



Second, the equation for T 0tl (conservation of energy -zeroth component of the 
four-momentum) takes the form 

-4 / vVc) 2 + |q| 2 / dq + -dW x f qfdq = 0. (15) 

to at j K 3 m j R 3 

And third, the equations for T 4AI , which represents the conservation of the spatial 
part of the four-momentum, is 

-! / qVrfq +-#- / , dq = 0. (16) 

to dt 7 R3 to dxj J R 3 v?(mc) 2 + |q| 2 

2.3 The H-Theorem and the entropy 

An important property of the BGK model (fTU]l is the H- theorem, which gives 
the local dissipation law for the functional L 3 /ln(/) <iq, namely, the following 
equation: 

^d t f f\ n (f/ V )dq+^dw x f q/ln(//T,)^ < 0. (17) 

w J R 3 to 7 r3 g u 

In order to achieve the iJ-Theorem, it suffices to show the dissipation property 

- / < 0. (18) 

m JR3 q l 
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This is already proved in [TH1 Chapter 8, eq. (8.8)]. 

Let us see next that the local Jiittner equilibria arising in our framework 
are determined by an extremality principle that is akin to that of maximum 
entropy. 

Lemma 2 (Minimum free energy principle). Let /(q) > be given and let Jf be 
the associated Jiittner equilibrium constructed in Lemma\J] Then, the following 
inequality 

k B (3 \ ( k B P 



(mc) 2 J j \ (mc) 7 f 



e > 



holds. 



Proof. The proof is similar to that of Proposition 2.1 in [40]. Using a Taylor 
development, we get the identity 

Mf/v) = JMJf/v) + (f - J f )MJf/v) + 1] + (/ ~ 2g J/)2 , 

being £ some momentum distribution that can be written as a convex combina- 
tion of / and Jf. Then, taking in to account Lemma [TJ the difference between 
the entropies associated with Jf and / can be written as 



o-jj-c-f = ^(u/)|. 



q nn(J f /r,)(f-Jf)-±+ / q 



dq , f _„ (/- J» 2 rfq 



m \r,(mc) 3 M(Pf)J 1 H U /; q° 

--%(«/). / MV-JfA 



T 

Taking into account (JT]) we deduce that the last term above is non-negative. 
The first term is already zero thanks to Lemma [TJ Next we use that j3j f = (3f 
to write 

<rj t -<Tf>- — ^(Uf),(uf) u [ q» l f'(f-J f )*L=-** {fi p } _ 
Altogether the result follows. □ 



3 Asymptotic limits of the relativistic BGK model 

One of the aim of the paper is to perform three distinguished limits of the BGK 
system (j 10|) : the non relativistic limit, the hydrodynamical limit and the ultra 
relativistic limit. In order to do that we will consider the physical meaning of 
the involved constants represented on a set of only four typical variables, which 



allow to rewrite (fTU)) in a dimensionless form for its mathematical study. We 
will precise later the connections between these four variables which give rise to 
the three limits. Let us choose a scaling such that: 

q = q/i, x = xL, t = It, 

where "~" stands for dimensionless numbers and \x, L and r stand for the typical 
microscopic momentum, the typical length and the typical time, respectively. 
Also, we perform the change of unknown 

f(t,x,q) = ^f(I,x,q), 

where N is the typical density. Then, system (fTU)) becomes now 
d ^ l ru\ q Vs/ , , 1 



j ( Zh_ \ q Vs ^ 

dI J + \mL) I |JL ,2 



(cur) 



(19) 



I mc 1 

where = ^ «// must be specified for each limit. Also, (|17p becomes 



ft/ /V^W 



( — F )div 5 



q/ln( 



< 0. (20) 



i+ — q 

I mc *| 

Then, the three limits will be attained by choosing the following relations: 



Limit 


Scaling relations 


Asymptotics 


Non relativistic 


mL 
M= 

T 


LOT = CSt 




fi « mc 


Hydrodynamic 


mL 
M= 

T 


L 

c = — 

T 


77/i 3 

AT =cst 


1 

— << r 


Ultra relativistic 


LOT = CSt 


L 

c = — 

T 




TOC « jj, 



(21) 



Also, a dimensionless relativistic BGK model will be analyzed in Sections |4] 
and [5] when the nondimensional constants involved in (fT9")) -(f2"U )) are 1, which 
corresponds with the choice 



Normalized scaling relations 


L 

c = — 

T 


LOT = 1 


N 


/i = mc 



(22) 



3.1 The non relativistic limit 

In order to give consistency to our approach to the relativistic BGK model, 
let us recover in the non-relativistic limit the classical BGK system, see also 
[17j[20j[59] for a similar analysis in this context. 
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In order to perform the non relativistic limit we choose the first scaling 
given in (f2Tj) . It corresponds to the assumption that the typical microscopic 
momentum fi is small compared with mc, i.e., the speed of particles is slow 
compared with the speed of light. Then, setting Q := lot and fj := r]/j, 3 /Af 
-assumed to be constants- and defining e := /i/(mc) << 1, our system ([T! 
becomes now 



difn 



1 



sA + l £ ql 2 



Q ' V 7 xfnr — 



(23) 



where it only remains to write in a compact form the expression of the right 
hand side. To do that wc write: 



1 



1 



l-A + O^)), 



(24) 



and then, after the change of variables dq = fi^dq on the integrals, we get 



n j 



[ /(i,x,q)dq) - V" ( / cqj(i, x, q)^§ 

/ fnr{i,x,q)dq) -e 2 ^2( qif nr (i,x,q)dq 



0(e 4 *25) 



Defining then the classical (scaled) density as usual, 



/„,.(*, x,q)dq, 



we obtain the asymptotic behavior of n/ 

rif(t,x) = Nn nr + 0(e). 
In a similar way, we use (|24|) to find 

dq W> f qf nr (i,x,q) 



(26) 



cq/^.x.q)-^ = 

q v mc 



V^+7 2 H 2 



dq = eAf qf nr dq+0(e 2 ) 



which in combination with (|26|) allows to get 

— / qf nr (i, x, q)dq) + 0(e 2 ) := eu„ r + 0(e 2 ), 



(27) 



where u nr is actually the (scaled) classical velocity. In order to obtain the 
expansion for (3f we first combine a Taylor expansion of — s 2 around s = 
together with ([25]) to obtain 



- = T7^(l + ||u n ,.| 2 )+G(e 4 ). 



(28) 
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Second, we use to obtain 

/ f(t,X,Cl)^ = ^( [ f nr (i,X,q)dq-e 2 [ ^-/„ r (i, x, q)dq) +G(e 4 ), 
Jr3 g° mc 2 \Jr 3 Jr 3 2 / 

which, combined with (|28|) and using the definition (|13[) of j3f, allow us to write 

|^(/3/) = l + f (W| 2 - — f |q| 2 /„ r (i,i, q )d q )+0(e 4 ). 
&2 2 \ n nr J R3 J 

In particular, we observe that K\{[i j) / K.2{(3 j) approaches 1, so equivalently /3/ 
is large. This is precisely why we use the following expansion (check |19|). valid 
for pf » 1 



Finally, we use the two expansions of Ki/K 2 to obtain 

^- = i(— ( |q| 2 /nr(i>,q)dq- |u„ r | 2 ) +0 

Pf 3 \n nr J R3 J \ Pf . 

and, actually, define the (scaled) temperature at equilibrium, 



T nr = ^--=\(— I |q| 2 ./W(i"S,q)dq- |u„ r | 2 V (29) 

Pnr o \n nr Jjfs J 

which corresponds with the classical one. With this notation we can also write 
the expansion for pf 

e 2 p f =P nr (l-(D(e 2 )). (30) 

Let us now check the behavior of the Jiittner function in the non-relativistic 
limit. We first observe that (see pQ) 



1 ePt 



Mifif) 



3/2 



H(£))-£(£rV«*»- *> 



On the other hand, using (|27|) . we can get the following equality for the size of 
the exponent in the Jiittner function 

= ^l + | U/ |Vl + c 2 |q| 2 - eu, • q = 1 + ||q - u Ilr | 2 + 0(e 3 ). (32) 

Then, putting all the expansions {25]), {27]), (|30j!, ([3TJ and (J32J together, we 
finally obtain 

Jf Jf = ^ (IF) 3/2 ° XP {- ^~2 nr? } + ° {e2) := G{fnr) + ° {e% 
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which is, up to order e 2 , the classical Maxwellian G(f nr ), that can be also 
written as 

Then, (|23|) becomes (up to order e 2 ) the classical BGK equation 

dtfnr + q • V x / nr = Uj(G{f nr ) — f nr ). 

We can conclude that, assuming the convergence of all the involved moments of 
fnr, the limit satisfies the classical BGK equation. 

Next, for the energies, identities (j2"l)) and (|52"j) lead us to the following ex- 
pression 

e f = Nc 2 (n nr (l - ^-|u nr | 2 ) + e — J ^ |q-u nr | 2 /„ r (t,x,q)dq^ +C(e 3 ). 

The previous relationship can be seen as the usual expansion of the relativistic 

expression for the kinetic energy for small velocities: the first term, J\fc 2 n nr , is 

the rest mass energy of the system, while up to first of order in e 2 , ^-n nr \u nr \ 2 

is the (Newtonian) kinetic energy of the center of mass and \ J R3 |q— u rar | 2 f nr (i, x,q)dq 

the internal (Newtonian) kinetic energy of the individual particles (with respect 

to the center of mass frame). The latter is a measure for the temperature of the 

system, such that with (|29|) the above expression can be written as 

Afc 2 / 3 \ 

e/ + e 2 — n nr \u nr \ 2 = Mc 2 [n nr + -e 2 n nr T nr j + 0(e 3 ). (33) 

Note that in the rest frame (u nr = 0) this expression reduces to the well-known 
relation between the kinetic energy and the temperature of a monoatomic ideal 
gas (except for the rest mass energy term). 

Finally, we can also check the relation between the relativistic and classical 
entropies. Using again the previous expressions for Uf, we find, 

= _knM I fnA ^ )ln r inr(t,_x,q) \ d q + 0(£2)j (34) 
™ Jw> \ V ) 

where the first term of the right hand side corresponds, up to order e 2 , to the 
classical entropy. We also stress that, in the light of what was done in here, 
Lemma [2] reduces in the non relativistic limit to the usual maximum entropy 
principle for the Maxwell-Boltzmann distribution. This is a direct consequence 
of relation (|34| and the fact that the local energy reduces to the rest mass in 
the non relativistic limit -see (|3"3")) -. but then both the density of / and that of 
the associated Maxwellian are equal. Thus (fie)j, — (/3e)/ cancels in that limit 
and the functional in Lemma [2] converges to the classical entropy. 

3.2 Derivation of the relativistic Euler equations 

This part is devoted to the connection between the relativistic BGK model 
pop and the relativistic Euler equations. These are obtained as a particular 
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limit of the BGK model. Obtaining hydrodynamical hyperbolic limits from the 
relativistic BGK model is one of the goals of this paper, but the analysis given 
in the previous subsection allows to connect this type of limits in the relativistic 
setting with their counterpart limits for the classical BGK system. At least in 
the framework of regular solutions, hyperbolic and non relativistic (Newtonian) 
limits should commute. Then, in this sense we can assure that the relativistic 
Euler equations obtained in the next Section [3. 2 1 are the relativistic counterpart 
to the classical ones derived from the classical BGK model. 

In order to proceed with the hyperbolic hydrodynamic limit, we choose the 
second scaling given in (|2"Tj) . It corresponds to a fluid behavior: we assume 
that the time between collisions 1/lo is small compared with the typical time 
r, and the fluid description takes place. Then, if we take now e = lot, we 
put r]/! 3 /AT — 1 in order to abbreviate the expression for the entropy and we 
eliminate the superscript "~" to simplify the notation, we obtain the following 
fluid-scaled version of the relativistic BGK model: 

dtfhu + q ■ Vx/k = —^—{J fiu - fhJ (35) 



where now q = q/ \Jl + q 2 and Jp stands in dimensionless form: 

J f^ = JT§iD" v (v /l + u ^ v/TW - u/ - ■ q ) } 

Also, the conservation laws (TT4")) , (fTS")) and (fTB")) read now 

9 f « .,_ , J: _. f q/k 



-/ f Eu dq+div x ^Ulf= dq = 0, (36) 

ot Jr3 J r3 y/l + |q|2 

f y/1 + |q| 2 /k^q + div* / qfhudci - 0, (37) 

|/ qV^dq+A/ ^^ dq = . (38) 
ot Jr3 dxj J R3 y/i + | q |2 

Let us denote by {f% u (t, x, q)}e>o a sequence of nonncgativc solutions of the 
equation (|35| such that f Eu converges almost everywhere to a nonnegative func- 
tion Jeu, as e goes to zero. Assume also that all the moments converge to their 
corresponding moments and that Jf^ converges to Jf Bu in an appropriate sense 
to be specified later. Then, multiplying (f3"5|) by e and letting e go to 0, one has 
in the limit: 

fEu = Jf Eu - 

Then, deriving macroscopic equations from (|36|) . (|37j) and (|38|) . will require 
information on different moments of the equilibrium Jf Eu , stated in Lemma [121 
Under suitable convergence assumptions (which we prove below), we can ensure 
that 

f A(q)f Eu dq^ f A(q)f Eu dq= f A(q)J fEu dq, 

Jr 3 Jr 3 Jr 3 
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where A(q) stands for all the involved moments of Je u , that is, 



A(q) = 1, q\ y/T+W, 



q 4 q'q^ 



y/T+W y/T+W 



These convergences combined with (fBl)]). ([37]), ([35]) and the computations of 
Lemma [12] lead to 



ft [nf Eu \fl + |u /E J 2 ^ + div x (n fEu u /b „ ) = 0, (39) 

ft (-i>/ B „+^/ B «x(AfaJ(i + l u / B J 2 )) 

+div, (n fEuX {Pf E Juf E „y/l + |u/ B J 2 ) = 0, (40) 



|u/ E J 2 



+ft, (p/^ + ^„x(feJu} E „u; E J = 0, (41) 

respectively, where xW) — j$ + ^(Z?) with \&(/3) given by (|59[) in the Appendix. 
The result that we have sketched can be stated as follows: 

Theorem 1. Let f Eu (t,x,q) be a sequence of solutions of the equation (|35[) 
with nonnegative initial condition < /J; u (0, x, q) = /^(x, q) verifying that, 

( ( (1 + VTTIq7+N + |ln(/^( a; ,q))|)/^( a ;,q)dqd a; <C<oo. 

Assume also that f Eu (t,x,q) converges almost everywhere to x, q), as 

e goes to zero. Then, the pointwise limit 2, q) is a Jiittner distribution 

whose moments are the corresponding limits of those of f Eu . In particular, the 
functions nf Bu ,Uf Eu and fif Eu associated with fEu solve the relativistic Euler 
equations (|39[) (|4ip . Moreover, the limit fE u {t,x,c\) does also satisfy the fol- 
lowing entropy inequality 



ft (n fEu Jl + |u feu p (K M "^ ) ) ~ ft».*G3/*J 

+div, (n /Eu u /E „ (ln( M ^ < <?42) 

Proof. We first prove that the solutions are also nonnegative functions. To 
do that it suffices to observe that the Jiittner function on the rhs of ([35]) is 
nonnegative and then 

^ /!;«(*, a + 9*, q) >-/!«(*, a + gt, q) 

and so, integrating, we deduce f% u {t, x, q) > e~ t f E l(x — q) > 0. Second, as 
in the classical kinetic theory, we will prove the following a priori estimates 

/ f {l + ^l+W+\^\ + \Hfhu(t,x,cO)\)fEi(^<l)dc l dx<C(K)<^, 

Jr 3 Jr 3 
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valid for any t in an arbitrary compact set [0, K]. The first two estimates on 

/ / fhu( t ^,q)dqdx and / / y/TT\kTfEu( t > x >$ < fa dx 

Jr 3 Jr 3 Jr 3 Jr 3 

follow directly, using the non-negativeness, from and (|3T|) respectively after 
integration with respect to x. Also, multiplying (|36[) by |a;| and integrating, we 
obtain 

4- f f Mf Eu dxdq= [ f f Eu dxdg<C. 

dt J R 3 J R 3 J K 3 J ffi3 |xk/l + Iql 2 



Then. 



/ / \x\f Ea (t,x,q) dxdq< f f \x\f e / u (x,q) dxdq + KC, VtG [0,if]. 

JR 3 JR 3 JR 3 JR 3 



On the other hand, (|20|) reads now 

ft/ fhMfhJ dq + d\v x f — g— /IJn^Jdq < 0, (43) 

JR 3 JR 3 V 1 + l^l 

and so, integration on i 6 R 3 gives us 

/ / /kM/k) <fa rfq < / / fEuHf E i) dxdqVte [o, if]. 

JR 3 JR 3 JR 3 JR 3 

The last a priori estimate on /J;,Jhi(/J; u )| follows straightforwardly from the 
following trick due to Carleman: 

s|ln(s)| < sln(s) + ks + 2e~i, Vs > 0, Vfc > 0. 

by choosing s = f Eu (t,x,q), k = \x\ + y/l + |q| 2 and using previous estimates. 

These a priori estimates allow us, via Dunford-Pettis theorem, to pass to 
the limit weakly in L 1 (R 3 x K 3 ), uniformly in t G [0, if]. The pointwise conver- 
gence together with this weak convergence allow to pass to the limit strongly 
in i 1 (K 3 x R 3 ), uniformly in t G [0, If]. Then, the associated quantities n/j, , 
u/« and fin pass to the limit strongly to their respective ones. As J/ Bu is a 
continuous function on the variables n, u and fi we get that Jp — > Jf Eu ■ On 
the other hand, as we said before, multiplication of (|35j) by e and passage to the 
limit for e — > give that the limits of J/« and that of f Eu coincide, at least in 
a distributional sense, then Jf Eu = Jeu a.e. 

Finishing the proof only requires to take limits in ([3"6" |) - (j3"5]l to obtain (I3U1) 
(|4Tj) and in (l43l) to obtain (|42|) . The only two involved quantities that pass only 
weakly to the limit are 

/ q/k(* 5 a,q)dq and I f Eu (t,x,q)ln(f Eu (t,x,q))dq, 
Jr Jr 
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but this is enough for our purposes, because the pointwise convergence allows 
to identify their weak limits as 

/ qJ/ E „ (*, x, q)dq and / J fsu (t, x, q)ln( J /e „ (t, x, q))dq, 
Jm. Jm. 

respectively. So, computations of Lemma [TBI concludes the proof. □ 

Remark 1. As in the classical case, the connection between kinetic and macros- 
copic Euler fluid dynamics results from the properties of the collision operator. 
These are the conservation equations and the entropy relation (which implies 
that the equilibrium is a Maxwellian distribution for the zeroth-order limit) 
[21 03 E2f- Analogously in the relativistic case \40[ \57j these properties allow 
to obtain the relativistic Euler equations as an hydrodynamic limit. Therefore 
it is very natural to obtain here the same relativistic Euler equations. 

Remark 2. Let us mention that the equation of state that we would obtain in 
the formal hydrodynamical limit to the relativistic Euler equations -a state equa- 
tion which consists on extrapolating the relations among the thermodynamical 
quantities of local Jiittner equilibria to arbitrary solutions of the fluid equations- 
is already discussed in \57^j , under the name of kinetic equation of state. The 
results in \57\ \18\/ show that under such equation of state the relativistic Euler 
system is hyperbolic and causal, and the speed of sound is bounded by c/y/3. 



3.3 The ultra-relativistic Limit of the relativistic BGK 
model 

Let us now consider the ultra-relativistic limit of our relativistic BGK model. 
Some of the results in this sense -as the limiting behavior of the Jiittner equi- 
librium or the ultra-relativistic Euler equations- have been partially discussed 
elsewhere ( [El [39l 00] ) , but we include them to give a completely coherent pic- 
ture. The ultra-relativistic limit corresponds to the third scaling given in (|21j) . 
that is rac « fx, and can be motivated as an asymptotic behavior when the rest 
mass m is very small compared with /i/c, as it happens with neutrinos or even 
photons in the limit case of zero rest mass. Then, normalizing the remaining 
two constants ui := ujt and fj := nfi 3 /Af as in section 13. 1[ and defining now 
e := mc/fj, << 1, the system (fT9|) becomes 

dtfur + T^TW = V^TW ^ Jf " Ur ) ' m 

Here f ur = /i 3 f/Af denotes the scaled distribution under this ultra-relativistic 
scaling. As in subsection 13.11 it essentially remains to work on the expression 
of the right hand side. To do that, we define the ultra-relativistic quantities 
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associated to the distribution /„ 



f ur (t,x,q)dq) -Y"(/ -r^fur(t,x,q)dq 
J r-r \Jk 3 111 



u„ 



A, 



■rrT/ U r(i,x,q)dq /ji, 



T=rfur{t, x,q)dq /n v 
3 IQI 



(45) 



We argue analogously as in section IQI to obtain the following asymptotic ex- 
pansions 



1 



"I 



1 



1 



A*c ^e 2 + |q| 2 Mc M 

2 3 



AT 



A/' 



q 



qi/«r^q 



3 \/e 2 + |q| 2 y 

/ ur (F, x, q)dq = u„ r + £>(e 2 



(46) 



= Mn ur +0 (e 2 X47) 



(48) 



Also, to obtain the behavior of /?/ we first use (fl~3|) combined with the expansions 
(1461) and (H3 and definition (EH to deduce 



*T 2 



08/) 



dq 



1 



rif J R a y>e 2 + [qp n„ r 7 R 3 |q 
^/3 ur + 0(e 2 ). 



/ ur (<"x,q)dq + 0(e 2 ) 



In particular, note that now Ki(f3f)/K2(f3f) approaches 0, so equivalently j3f is 
small and justifies the following expansion, which is valid for /?/ << 1 (see [T]), 



/?/ 



g(/3 / ) = ^ + 0(/3 2 ). 

Using finally these two expansions of K\ /K 2 we conclude 

P f = e/3 ur + 0(/3j) + 0(e 2 ). (49) 

We are now ready to analyze the behavior of %-Jf- We first write it as 

ex P { - ^(^l + lujlV^TlqF-u/ • q) } 



Tr J f 



M 3 Afe 3 M(p f ) - i e 
and then, using the asymptotic relation (see [T]) 
1 1 



M(/3/) 8tt 



= — Pf + 0(!3}) for/3 / «l 
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and the expansions (|37]l. (gSJ) and P9")l , we deduce that 

^J/ = ^^exp{ - ^(iqlVl + luurP-tW • q)} + 0{e). 
Then, up to order e, (El) becomes the ultra-relativistic BGK equation 



$tfur I i — i ^ xfur 

|q| 

Q ( n ur /3: 



q| V 8tt 



expj - /3 ur (|q| >/l + |u ur | 2 - u ur • qj i - /„ 



To conclude we compute various macroscopic quantities of local equilibria in 
this limit. The calculations are based on the asymptotic expansion of \&(/3) (|59[) 

*(/3) ~ | + 0Q3) for/3«l (50) 

(see [T] for instance). Thanks to Lemma [TT1 we have 

2 r/fl , 3 2 ^ 
ej. = c ri/W(p7) ~ — c n ur — := e ur — . 

' Pur e e 

In a similar way, 

o Tl f 2 n ur A/" A/" 
Pj/ = c ~ c := p„ r — . 

P/ Pur £ e 

Note that in this regime the energy amounts to three times the pressure. We 
can also use Lemma [TU] to obtain 

N N 



~ Nn ur < r , T»» p ur g^ + A— Pur < r u ur 

1 6 6 



4 Linearization of the relativistic BGK operator 

From now on, we will use the dimcnsionlcss system (|19[) by using the choice 
|). If we skip the superscript "~" to simplify the notation, it becomes 

dtf + q-V x f = ^(Jf-f), (51) 



where q° = \J\ + |q| 2 , q = q/q° and the local equilibria given by Jiittner 
distributions (or relativistic Maxwellian) takes the dimensionless form: 



J(n, p, u; q) = " exp{-/?(V(l + |u|*)(l + |q| 2 ) - u • q)}. 
We will study solutions close to a family of global equilibria, those of the form 
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From the mathematical point of view it is handy to bring in a new variable a 
defined by means of 



] (a), where a := — [ 



It is in one-to-one correspondence with f3, see Section lA.41 Many related formu- 
lae can be written down in a more compact way using this parameter instead. 
Thus, we hereafter consider 

-X( ao )y/l + \q.\ 2 -X( ao )q° 

J ° s J ^ ° ; = M(X(a )) = Mfxwy ° <a ° = X " {M < 1 



Our aim will be to find an equation for / such that g = J° + \J®f is a 
solution to the relativistic BGK model (f5T|) . To do that we will linearize the 
relativistic BGK operator around J° by using Taylor's formula. It will be useful 
to control the moments of the perturbation. 



Remark 3. In the following we shall consider perturbations h = V J°f which 
need not be positive, so that objects like nh become meaningless. Nevertheless, 
we notice that in the following the object does only appear involved in com- 
binations like n\, nh,Uh, nhyfl + \uh\ 2 or n^ah- These are just integrals which 
make perfect sense even if h takes negative values. So, in order not to distinguish 
the cases in which h > from those in which this is not so, we define 



"1, 



»*)'-£0L*f)"- 



n h u h := I qh Q 



</q 

7 



n h y/l + \u h \ 2 := / h dq, 

it 3 

Jr3 q° 



and in this way we can keep the same notation as in the nonnegative case. Let 
us also remark that nhVLh, nh\J^ + |u^| 2 and nuau are linear in h. 



Lemma 3. Let g = J° + h with h = fyj®. Then, the associated macroscopic 
quantities to the function g are given by the following formulae: 



(*) ri g = J(l + n h ^/l + |u h | 2 ) 2 - n 2 h \u h \ 2 = ^\ + n 2 h + 2n hy /l 

(U) UgUg = U h U h , 



|u,| 2 , 



(in) n g a g = { = n h a h + a . 

jr3 q 
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2 

2 1 



Proof. It is mostly straightforward. We first note that 

f i Tti , s ^ /" , dq 
n g u g = / q(J = / qh— = n h u hl 

which proves (ii). Now we obtain 

n a g (l + \u g \ 2 ) = fj^(J° + h) = ( 1 + J R , hdq ) = (l + W! + I«hl 
which combined with (m) gives (i). To deal with (m) it suffices to write 

n 9 a 9 = / +ft)-^=Qo + ifia fl . 
Jm.3 Q 

This completes the proof. □ 

Let us denote by D 2 * J the 5x5 Hessian matrix of the Jiittner function 
with respect to the variables (n, u, a), at the point (ng, ug, ag) := 9(n g ,u. g , a g ) + 
(1 — 0)(1,O, ao), for any given 9 e [0, 1]. Being given J g = J(n g , a g , u g ; q), one 
has the following Taylor expansion. 

Lemma 4. Let g = J° + h be a solution of the relativistic BGK equation (|51[) 
with h = fy/J®. Then the following equality 

%_ 

n „ 



J 9 -J° = K-l)J (q) + ^u,-qJ°(q) 



+ / (1 - 6){n g - l,u g ,a g - a )D 2 J e (n g - l,u g ,a g - a )d9 
Jo 



Tj + T 2 + T 3 + r 



holds. 



Proof. The proof is a direct application of Taylor formula together with previous 
LemmaEl after computations of the first derivatives of the Jiittner function with 
respect to n, u and a, which read 

dJ J „ T a ( q° \ T 

-5- = -, V u J = /3 q u J, 

<9n n \^ ^1 + |u| 2 J 

respectively. Note that we only need to use the above formulae for the values 
(n, a, u)= (l,ao,0), which simplifies the expressions. □ 
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For the next step we set 

3a o . 2 i 
k := — h a - 1. 

Po 

Recall that g = J° + /Vj° = J° + h. Now we let P(f) be the following linear 
operator (see Remark [3]) 



ao<7° - 1\ /^— -; ,2 



K / \ K 



*(A>) - 9 



VWoK - v/i + l^l 2 )^ ./Jo + / »h(aoy/l + |u h | 2 -^) \ 



P(f) 



\ «0 / V K ° 

+/3 n h u h • q\/j° := A f VJo + B f q VJ° + C f ■ qy/J°. (52) 

The following result shows how does the operator P enter in our framework. 

Lemma 5. Let g = J° + h be a solution of the relativistic BGK equation (|5ip 
with h = fy/j®. Then f verifies 

d t f + q-V x f=±(P(f)-f)+T(f) 

where the nonlinear term F(/) is given by 



2 4(1 + 2V2 + y/l+V h )J q° 



a -Dh VJ° 

-PQ-, ■ — , nh\ih ■ q — — 



Vh VJ° 

' nh0h (l + VT+V- h )VT+V- h q° 

(nj { V h 3 -3V h 2 \ <f(Po)-q° ^ 

~ a ° V 2 2(2 + V h -V h 2 + 2y/l + V h ) J kq ' q° 
1 f 1 

- n r-^ / (1 - 0)(n g - 1, u s , a g - a )D 2 J (n g - 1, u g , a g - a )d9. 
q u V J u Jo 



Here J := J(9(n g , u s , a g ) + (1 — #)(l,0,ao)) and T>h is given by y/T+Vh 
n g = yjl + n 2 h + 2n h ^/l + \u h \ 2 . 

Proof. We first remark that if g = J° + f\f~J® verifies (|5"Tj) . then / satisfies 

1 ( J 9~ J ° 

and the proof consists in identifying the linear and the nonlinear parts on the 
right hand side. Actually, we have to prove that 

y(^-/) = yW)-/) + r(/). 



d t f + q-V x f=-^(^ T =^-f 
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In order to do that we use Lemma H] to decompose J g — J° into its linear and 
nonlinear parts. First, we use the following equality 



v h v h 2 



2 A(l + V h /2 + VT+Vh)' 
to decompose T\ as 

Ti = (y/l + V h -l) J° 



+ 2n hy /l + \u h \ 2 V h 



2 



4(1 + V h /2 + VT+Vh) 



J° 



v h 2 



V h /2 + VT+V^) 



J° := Lx+Ti. 



Second, by using the equality 

1 . V h 



(53) 



we can write T 2 as 

T 2 = (3 Q n h u h ■ q J° + /3 ° :P ' 1 ^ «/.Uf. ■ qJ° := L 2 + T 2 . 

\/ i + V h (i + ^ 1 + U h ) 

For the third term, we first note that M'(/3)/M(f3) = -*(/?), sec flSSJ) in the 
Appendix, and also the fact that 



so that Af' (ao) = 1/kq. Then, T3 can be rewritten as 



Using again (1531) , we can write T31 as 

g(go) - g° o Vt 

r 3 i = — n h a h j + /ttwtTT 7tTtt\ n ^ J '■= L 3 + r 3 . 

Finally, using the equality 

1 __X> h V h 3 -3V h 2 



VT+Vh" 2 2(2 + % - 2?/i + 2V1+X>h) ' 
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we can write T32 as 
T32 = a 



n\ + 2n h ^l+ | U/l | 2 



2(2 + V h - V h 2 + 2^/1 + V h ) J ko 



-a — n h y/i + \u h \ 2 J° 



Vh* - 3V h 2 \ g(ft) - g° 70 

2(2 + X> h -2? /i 2 + 2Vl+^/ t )7 «o 



We only have to remark that P(f) = (L\ + L2 + Lj, + L^j/vJ and that 
r (/) = ( r i + T 2 + r 3 + r 4 + f)/(q°y/j°). Just note that 

# (AO - 1° _ a oq° - 1 



1 - a - 



helps to deal with the factor of Uhy/1 + |u/j| 2 in the sum of Li and L 4 . □ 

4.1 Analysis of the linearized operator: the projector over 
the distinguished space 

In this paragraph we determine some properties of operator P that will be 
needed for future analysis. We could consider the interplay of the projector 
with different scalar products. At least two of them come quickly into mind: 



(/, 9) 



I /(q) S(q) rfq and (/, g) qo = f ZMlM dq . 



To proceed, let N be the five dimensional space given by 

N = span{VJ°,q IJ 'Vj°}. 

Lemma 6. The linear operator P given in (5ty) is the projection from L 2 (M. 3 ) 
onto N with respect to the scalar product (-, -) qo . In particular, it is self-adjoint 
with respect to that scalar product. 

Proof. We first observe, from the definition (|52")l of P(f), that it is a linear 
combination of the basis in N and that 



P(VJ°) = V JO, P(qWJ°) = qWJ°, P(q VJ°) = q V J°, 



which is true by using a direct computation and the fact that l—^^(fio) 



■ 

Then P acts on TV as the identity. 

On the other hand, we can easily compute, for any / and g, that 



(P(f),9) 



«o po + qo , *(/3o) , a p q 



k nopoqo nopoqo poqo 



A) — /(q)3(p)v^°(q)v / >(p)^P 
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Let us define Lf = { J, as the linearized relativistic BGK operator. We 



which is obviously a symmetric expression and shows that P is self-adjoint. 
Finally, we write for any / and any g £ N 

(f-P(f),g) qo = (f,g) qo - (P(f),g) qo = (f,P(g)) qo - (f,P(g)) qo = 0, 

which concludes the proof. □ 

P(f)-f 
have the following result 
Lemma 7. The operator L satisfies the following properties: 

1. It is self-adjoint with respect to the scalar product (•, •). 

2. Ker(L) = N. 

3. L is non positive. In fact, 

(Lf,f) = -((I-P)f,(I-P)f) qO <0. 

4- L is decomposed as L = — 1/q Id + K, where K is the compact operator 



in L Z (W) given by K = 

Proof. By definition, Lf = is equivalent to f — P(f), so clearly / £ R(P) = 
N. The sclf-adjointness of the linear operator L is easily obtained computing 
from the definition (1521 that 



■" /r3 \k K po K q" n p q v p <7 u 



JR3 q v 



To prove 3 we first notice that 

c I /./)-<en F Z.pcfl-/> + <^P.pcfl>-(^./>. 

Therefore, using Lemma |5] we observe that the last two terms cancel and this 
completes the proof. 

The decomposition of the linear operator L is trivial. Note that the oper- 
ator K = P/q Q is compact because its range R(P) = N has finite dimension. 
Moreover, for future development, let us write it in the Hilbert-Schmidt form 
as follows 

K{f) = / fc(q,qi)/(qi)dqi, 
where the kernel k is given by 



«(q,qi) = i i- 



qo iv k 
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Thus, K defines a Hilbcrt-Schmidt operator in L 2 (Rq), as it is easy to see that 
fc(q,qi)Gi 2 (i 3 xR 3 ). □ 



5 Existence of solutions to the linearized Rela- 
tivist ic BGK equation 

In this section we will show that the initial value problem for the linearized 
equation 

dtf + q- V,/ = Lf (54) 

has a unique weak solution in L 2 (R^ x R^) which is global in time. This is 
done by means of the semigroup representation of the solution, analogous to 
that given by [27j [50] for the non-relativistic Boltzmann equation. 

To proceed we need to introduce some notation. Let I > and let H l (x) 
denote the Sobolev space of L 2 (R^) functions, whose derivatives up to order 
I belong to L 2 (R^). We shall denote the partial Fourier transform in x of 
/ G L 2 (R 3 X x R3) as follows: 



/ (C, q) = (2tt) = / exp(-i(x)f(x,q)dx. 

JjSL 3 

Then, we set H l (() as the image under Fourier transform of the space H (x) 
with the following norm: 



ll/IUo= ( 1+ I C l)V(C) r ,„=\\f(x)\\ H 



I'(x) 



Finally, let Hi be the Hubert space L (W^H (x)) with the norm given by 



11/11/ = Jf M Jf(;<l)f H i {x) d q 



We define the operator B acting on Hi as 

B = L — . 

Being L a perturbation of the compact operator K, then its domain is given by 



0(S) = {/€fT«/^V./ + 



Now the equation (f5"4")) can be written down as ^ = Bf. Given / G Hi, we 
can take the Fourier transform with respect to the spatial variable in the above 
equation, thus obtaining 

9f _ a? 

m ~ Bf - 
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where 

C-q 



B = L- i- 



q° 



Following now the standard techniques for the non-relativistic case [491 EQ] , an 
important property of B is obtained. Namely, 

Theorem 2. For each (eK , the operator B generates a strongly continuous 
contraction semigroup on L 2 (Rq) such that, for any f E L 2 (Rq), one has 

exp(^)/(q)|| i2(K3) < ||/(q)|| i2(K 3 q) , f > 0. 

The operator B can be used to construct an explicit representation of the 
semigroup exp(£i?) for / E Hi. Using the same argument as in the non- 
relativistic case [27l [49] , one can show that the operator B generates a strongly 
continuous contraction semigroup on Hi. It is also shown that such semigroup 
is given explicitly as: 

exp(tB)f(x,q) = (27t)-l [ exp(i(x) exp(iB)/(C, q) d(. 

Moreover, the following estimate holds for any t > 0: 

||exp(tB)/(a:,q)|| J <||/(x > q)|| I . 
Then we can state finally the following result: 

Theorem 3. Let fo G H\ with I > 0. Then, there exists a unique f(t,x,q) 
global in time solution of J5^[ ) 7 satisfying f(t) € Hi and 

\\f(t)\\i < ll/olli, Vt>0. 



A Appendix 

The aim of this Appendix is to compute some of the various quantities involved 
in this paper in order to make it easier to follow. For the sake of simplicity, 
it is enough to perform these computations for the dimensionless quantities 
as in Sections H] and assuming the scaling ([2"2")l or, equivalently, to assume 
m = c = u) = r) = l. 



A.l Lorentz invariance 



Let A be a Lorentz boost (i.e. a linear isometry with respect to the Minkowsky 
metric) in R^. As we are dealing with particles of unit rest mass (the so-called 
mass shell condition) , this transformation A can be meaningfully seen as acting 
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on R q . Given any distribution function /, we can define a new distribution 
function f\ by means of 

f A (t,x,q) = f(t,x,Aq). 

As 

= (Av) f t(Az) t * for any v, z G K 4 , (55) 

we can check directly that n, e,p and a are Lorentz invariant, i.e. nj = nf A and 
so on, for any Lorentz boost A. Moreover, making use of the fact that the ratio 
^ is invariant under the action of A, we see that 

f q»f A (t,x,q)^ = f (A-Y)/(t,*,q)^J 
Jr 3 T Jr 3 Q 

= A- 1 (7 «?"/(*,*, q)^?Y 

This means that 

u /a = A_1 «/ 

(e.g. macroscopic boosts on the local velocity of the system are uniquely de- 
termined by the action of the same boosts -in a contravariant way- on the 
microscopic local velocities of the gas). The Lorentz invariance of the volume 
element ^§ shows also that j3 is Lorentz invariant, as the ratio 

Kitf) _ 1 f dq 
K 2 (0) n f J R3 g° 

is Lorentz invariant too. We summarize these facts as: 

Lemma 8. Given any distribution function f, the scalar quantities nf A , e/ A , 
Pf A , <jf A and (3f A are Lorentz invariant. The vector Uf transforms according to 
w/a = A ~ lu f- 

It is instructive to consider the special case in which the distribution function 
induces a local velocity that is found to be zero; that is, the physical objects 
that we are representing are at rest with respect to the reference frame that 
we use to describe them. This situation corresponds to distributions / having 
Uf = (1,0,0,0) -Lorentz rest frame. It is useful to display formulas for the 
macroscopic quantities of the gas in this the computations are simpler 

than in the general case and the results can be related to a generic distribution 
by means of Lorentz boosts. These read now: 

n f = N°= [ f dq, 

JR 3 

e f = T 00 = / VT+Wf dq, (56) 

JR 3 

p f = \\T^+T™ + T™] = \( |q| 2 /^. (57) 

•J O J R 3 <T 
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A. 2 Computation of the moments of the Jiittner equilib- 
rium 

We will need a more precise information about the moments of the relativistic 
Maxwellian. First we list for convenience some of them that can be easily com- 
puted in the Lorentz rest frame. Notice that in this case the Jiittner equilibrium 

reduces to 

J(n,/3,0;q) = " exp^^lTM^}. 

For future reference wc point that, using modified Bcsscl functions for the non- 
negative integer number j 

POO 

= / cosh(jr)exp{— j3cosh(r)}dr, 
Jo 

we can simplify some of the related formulae. For instance, we can write the 
function M (j3) given in ([5]) as 

M{p) = jKtf). (58) 
To simplify the notation we will introduce the function \t defined as follows 

Lemma 9. Let J = J(n, /3, 0;q). Then, the following equalities are verified: 
1. J d(\ = n, 



2 



. j q i Jdc l = f q i J^=0, 
Jm.3 J R 3 q" 



* L ^ £ - y> 

5. [ v/l + |q| 2 «7 dq = n^>{[5). 

Proof. The first relation follows from the very definition of M(/3) , and the second 
one just by a symmetry argument. To obtain the third one, we use integration 
by parts: 



H J ~K =47T / Jdr _ 

Ann [°° d / - 8 yT+p\„3 dr= [°° r 2 e - 0VT +^ dr 



e 



/3M(/3) J dr V J pM{fl) Jo 

3n f o /, t i„i2 , 3n 



/ e -/3Vi+lql 2 do 
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The fourth relation is a consequence of the following identity 

■dq = 47r/ ^ 2 dr=^-K x {p) 



v/l + |q| 2 Jo /? 

The sum of the third and fourth relations yields the fifth one by using (|59|) . □ 

The moments of the Jiittncr distribution in general form can be obtained 
thanks to the following classical decomposition (see [42] for instance): 

Lemma 10. The energy momentum tensor T^" can be expressed as: 

T^= [ q ^f^=-p f g^ + ( ef +p f )u^u v . 
Jm. 3 Q 



Then all the moments of a given / that appear as components of T Ml/ can 
be computed once we have the values of p/ and e/. This can be combined 
with Lemma [H which ensures that it suffices to compute the local energy and 
pressure in the Lorentz rest frame. These two are given by formulae (|56j) and 

(5ZJ. 

For the special case of / = J we can go further as the computations in 
formulae (|55|) and (|57p were already carried in Lemma HO Then we get the 
following result. 

Lemma 11. The quantities ej and pj are given by 

n 

Using the standard physical units, 

ej = c 2 n^((3), pj = c 2 ^. 

A direct application of the program sketched above yields then 
Lemma 12. Given any Jiittner distribution J , the following relations hold true: 

1. I q»J^=nu», 

2. ( |q| 2 J^= e , 7 |u| 2 + W (3+|u| 2 )=nvI/(/3)|u| 2 + (3+|u| 2 )^ 



v/TTIqlV dq = pj\u\ 2 + ej(l + |u| 2 ) = -|u| 2 + n*(/3)(l + |u| 2 ) 
4. q'Jdq = (ej +p J )^Y+Wu i = (n*((3) + £j y/T+^v*, 
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Proof. The first point follows from the definition of u M in terms of N^ 1 . For the 
remaining ones, we just take recourse on Lemma HU1 From there we get that 



f q°Jdq 

Jr 3 


= r 


[ q\J dq 

Jr 3 


= T 


/R3 <T 


= T 



'■T aa + T 33 =3p J + (e J +p J )|u J | 2 . 
This is to be combined with Lemma [TT1 To conclude, we notice that 

Jr 3 T Jr 3 t 
and the last relation follows. □ 

A. 3 Entropy fluxes 

We can compute also the entropy densities and fluxes of the Jiittner equilibrium, 
which are used to obtain information in the hydrodynamical limit, thanks to 
the H-thcorem. 

Lemma 13. The following equalities are verified: 

f ^^Jdq = u\pj + (ej+pj)\u\ 2 ), (60) 

Jr 3 q 

J^jHJ)d q = ny/rr^L(-^r\- ^{fs)\ (6i) 
/,i Jln(J)dq = mi H-m)-^)- 

Proof. Using Lemma [TU1 we get 

-^-Jdq = Uj(pjS ij + (ej + p J )u i u j ) = u\pj + (ej +pj)\u\ 2 ). 



/K 3 9 
To prove (jUTj) . note that 



J^J\n(J)dq = In^j^ j^J d q -(3^TTW J R3 V^+Wjdq 

+ flu- qj dq 

Jr 3 

and then using Lemma Q21 items 1 , 3 and 4 we obtain ([6"Tjl . 
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In the same way 

f ^j]n(j)dq = kx(-£-) [ ^jdq-p^T+W I q^q 

Jr3 q° \M{p)J J w3 q° J R 3 

+ P / -rrq • uJ rf q> 

and using Lemmas [11] and I12| items 1 and ^, combined with (|60|) we arrive to 
the last identity. This proves the Lemma. □ 

A. 4 Monotonicity of K 1 /K 2 

To begin with, let us recall the following recurrence relation: 

K 2 {P) = ^K 1 {p)+K {fi). (62) 



This can be used to show that 



KiW) ~ (3 
as K ((3) < Ki(/3). We also note that 



^ a C8)y /5K 2 (/5) + U2(«/ 



(63) 



3 /3 2 2 -/3 

> -„ — r 1 = 



(5 + 2 03 + 2) 2 (/3 + 2) 2 



which is strictly positive for /3 < 2. 

Next we analyze the case (3 > 2. For that we deal with the integral represen- 
tations of the incomplete Bessel functions. Using the substitution x = sinh(s/2) 
we get 



Jo \/l + x 2 



oo 



K (fi)+K 1 {p)= l~ (l + cosh(s))e- f3cosh ^ ds = 2 I 2 + 2x2 e -M+*?) dx 



By means of the inequality 



1 x 2 

> 1 for x > 



VTT^2 - 2 

we obtain the estimate 



V2^ / 1 



8^ 32^ 5 



1 We use that the binomial series is alternate and the fact that when we truncate the series, 
the error term has the same sign as the first term that is discarded. 
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Arguing in a similar way but using this time the inequality 



we arrive to 



Therefore 



128/3 



and 



Kpj^+K^) 256/3 2 + 32,3 - 24 
if (/8) ~ 128/3 2 - 16/3 + 9 



ifi(/3) 128/3 2 + 48/3 - 33 



if (/3) ~ 128/3 2 - 16/3 + 9 ' 

Notice that both the numerator and the denominator arc positive for /3 > 2. 
This estimate can be used in combination with (|52^1 to get 

K 2 (P) 128^ 3 + 240/3 2 + 105/3 - 66 
Kxifi) ~ 128/3 3 + 48/3 2 - 33/3 

We plug this into (|63|) so that 



#i(/3)Y >3 128^ 2 + 48/3 - 33 (128/3 2 + 48/3 - 33) 2 /3 2 



ivT 2 (/3)/ " 128/3 3 + 240/3 2 + 105/3 - 66 (128^ 3 + 240/3 2 + 105/3 - 66) 2 

_ 3(6656/3 4 + 8512^ 3 - 4080^ 2 - 2013,3 + 726) 
" (128/3 3 + 240^ 2 + 105/3 - 66) 2 ' 

Using that j3 > 2 we conclude with 

#i s ) Y > 3 ( 6656 /3 4 + 2419^ 3 + 726) Q 



if 2 (/9) / ~ (128/3 3 + 240^ 2 + 105/3 - 66) 2 
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